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Pixie Dust: The Silicate Features in the Diffuse Interstellar 

Medium 

J.E. Chiar^ and A.G.G.M. Tielens^ 
ABSTRACT 

We have analyzed the 9.7 and "18" fim interstellar silicate absorption features 
along the line of sight toward four heavily extincted galactic WC-type Wolf-Rayet 
(WR) stars. We construct two interstellar extinction curves from 1.25 to 25 fim 
using near-IR extinction measurements from the literature along with the sili- 
cate profiles of WR 98a (representing the local ISM) and GCS 3 (representing 
the Galactic Center). We have investigated the mineralogy of the interstellar sil- 
icates by comparing extinction profiles for amorphous silicates with olivine and 
pyroxene stochiometry to the 9.7 and "18" /im absorption features in the WR 
98a spectrum. In this analysis, we have considered solid and porous spheres and 
a continuous distribution of ellipsoids. While it is not possible to simultane- 
ously provide a perfect match to both profiles, we find the best match requires 
a mixture of these two types of compounds. We also consider iron oxides, alu- 
minosilicates and silicate carbide (SiC) as grain components. Iron oxides cannot 
be accommodated in the observed spectrum, while the amount of Si in SiC is 
limited to < 4%. Finally, we discuss the cosmic elemental abundance constraints 
on the silicate mineralogy, grain shape and porosity. 

Subject headings: ISM: lines and bands — ISM: abundances — ISM: dust — 
ISM: extinction — ISM: molecules — infrared: ISM 



1. Introduction 

Astronomical observations are often hampered by the confusing presence of foreground 
dust. These data are then adjusted using vaguely understood correction procedures, which 
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are perhaps best described as "faith, trust, and pixie dust." In this, astronomers follow 
notions which are widely spread in society because the term "pixie dust" is often used 
to refer to a technology that seemingly does the impossible. In this paper, we address one 
aspect of interstellar pixie dust: the nature of interstellar silicates and the infrared extinction 
"correction" procedure. 

The presence of silicate dust in the interstellar medium has been known for some four 
decades, but its precise mineralogy remains uncertain. Silicates in the interstellar medium 
(ISM) are thought to be the product of mass-losing oxygen-rich Asymptotic Giant Branch 
(AGB) stars. Amorphous silicates arc found in the circumstellar environments of almost all 
evolved oxygen-rich mass-losing stars (Kemper 2002). Crystalline silicates are found around 
AGB stars with high mass-loss rates, although it's possible these sihcates exist but are un- 
detectable around low mass- loss AGB stars (Kemper 2002). In the ISM, crystalline silicates 
are not detected and, recently, a stringent hmit of < 1% of the total sihcate population by 
mass has been placed (Kemper, Vriend, & Tielens 2004, 2005). 

Amorphous sihcates in the ISM exhibit absorption features with smooth profiles at 9.7 
and ~ 18 /im. These profiles can be fitted with extinction features based on laboratory 
derived optical constants. This procedure results in estimates for abundances in solid form 

of Si, Mg and Fe, which depend on the assumed porosity, particle shape, and the adopted 
silicate material (olivine or pyroxene-glass, Mathis 1998). Typically, such estimated require 
most or all of the available solar abundance of Si, Mg, and Fe. 

Bright late-type WC Wolf Rayet stars and IR sources at the Galactic center have been 

used to probe interstellar silicates in the diffuse interstellar medium (Aitkcn et al. 1980; Roche 
& Aitkcn 1985). The Galactic center lines of sight have the advantage of copious extinction 
{Ay ~ 30 magnitudes), but the disadvantage of being a crowded field where emission and 
absorption components can complicate the analysis of interstellar features. Late-type carbon- 
class (WC 8-9) Wolf- Rayet stars that are extincted by 6 to 12 magnitudes of interstellar visual 
extinction, provide an opportunity to study interstellar absorption features within the solar 
circle. Although these stars are known to have circumstellar dust, very httle of the visual 
extinction arises in the circumstellar environment (Roche & Aitken 1984). 

In this paper, we present the 2.38 to 40 /xm spectra from the Infrared Space Obser- 
vatory's (ISO) Short Wavelength Spectrometer (SWS) of four WG-type Wolf-Rayet stars. 
The data reduction is discussed in §2. In §3 we describe our method of fitting the spec- 
tral energy distribution of the WR stars. This estimated continuum is used to extract the 
silicate absorption features. In §4 we outline the calculations for the extinction profiles. 
These calculated profiles are compared to the best S/N observed profile in §5. Section 5 also 
compares our observed diffuse ISM profiles to observed silicate features in different environ- 
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ments. Section 6 gives an overview of the evolution of silicates from their production sites 
to their incorporation into young star environments. In §7 we discuss abundance constraints 
for the silicate minerals considered and in §8 we discuss other possible contributors to the 
absorption. Extinction curves from 1.25 to 25 /xm are constructed for the "local" interstellar 
medium and the Galactic center in §9. Finally, we summarize our results in §10. 

2. Observations and Data Reduction 

Wc present the 2.38-40 j^m spectra of the four latc-typc WC stars, WR 112, WR 118, 
WR 104 and WR 98a from the Short- Wavelength Spectrometer (SWS, de Graauw et al. 1996) 
on the Infrared Space Observatory (ISO, Kessler etal. 1996). These spectra were acquired 
with the Astronomical Observing Template (AOT) 01 speed 2 with an average resolving 
power (A/AA) between 250 and 600, and were first presented by van der Hucht etal. (1996). 
ISO-TDT numbers are given in Table 1. 

Off-line processing (OLP) version 10.1 data were reduced starting at the Standard Pro- 
cessed Data (SPD) level using lA"'. For each band, the twelve individual detector scans were 
normalized ("fiat-fielded") to a constant or a first order polynomial. Points that deviated 
more than 3 sigma from the average were cut, and the resulting spectra were then rebinned 
to the average resolving power for each band (Leech et al. 2002). Band 2C (7-12//m) spectra 
suffer from residuals caused by incomplete removal of the relative spectral response function. 
For these data, an extra reduction step using the "fringes" routine in lA^ was carried out 
after rebinning the individual detector scans. The WC spectra have low flux at long wave- 
lengths, so band 4 data (29-45 /im) have low signal-to-noise and have not been included in 
the final spectra shown in Fig. 1. 

3. Analysis 

The spectral energy distributions of the late-type WC stars are complicated by emission 
from heated carbon dust arising in their circumstellar environment (e.g., Williams, van der 
Hucht & The 1987; van der Hucht et al. 1996). Modeling their spectral energy distribution in 
detail is beyond the scope of this paper. Instead, wc will adopt a physical, albeit simplistic, 
model for their continuum dust emission in order to extract the interstellar silicate absorption 
features. 

Following Williams et al. (1987), we assume that the WR stars' mass loss can be modeled 
by optically thin spherical shells of amorphous carbon dust. We assume that the heating 
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of the dust is dominated by visual/UV photons from the central star and that the dust 
absorbs these with unit efficiency. The IR emission efficiency, Q (z/), of carbon dust is well- 
represented by a /9-law with {3=1 (e.g., Q{X) = Qo{i^/^oY'i Bussoletti etal. 1987, Mennella 
etal. 1995). The temperature is then described by. 



where To is the temperature at the inner radius Tq- The contribution to the flux by a shell 
is given by 

For the grain density distribution, we will adopt a power-law density distribution, n (r) = 
no{r / ro)""' ■ The flux is then given by, 

, (4 + /3) n,ro47rVg,i/^+/^-^ ^, 

= ,„ 3 „ ^{yo,yi) (3) 

i/n 



with 7 = (3 — a) (4 + /9) /2 and |/ = hv/kT^. T is the integral, 

where the integral limits apply to the inner and outer boundary of the shell. For a constant 
outflow a = 2 and j5 = 1, we have 7 = 5/2 When y,, is very small and yi is very large, 
the integral approaches the Rieman zeta function and is independent of the exact boundary 
values. The frequency dependence of the flux is then contained within the pre-integral factor, 
3 + /3-7 = 3/2. 

There are four parameters in the fit, a, /?, Tq, and Ti. We have adopted a = 2 and 
i3 — \ and allowed and Ti to vary, then evaluated the integral. We have used this model 
to "fit" continua to the observed spectra. The values of the parameter fits for each of the 
WR star spectra are given in Table 2. The model "fits" are shown in Fig. 1 and the resulting 
optical depth spectra are shown in Fig. 2. The optical depths, their ratios and the visual 
extinctions (from van der Hucht 2001) are listed in Table 3. The lower panel in Fig. 2 shows 
the spectra normalized to one at their peak in order to compare the 9.7 /xm profiles. The 9.7 
^m profiles arc similar from 8 to 11.3 ;um, while there is some variation in the absorption 
between features, from about 11.3 to 15.5 ;um, with the spectra of WR 104 and WR 118 
being the outliers. The central wavelength of the 18 \xm feature and the ratio of the features 
is similar for 3 out of 4 sources. The outher, in terms of the ratio of the silicate features, is 
the spectrum of WR 104. This could be due to uncertainties in the continuum that are a 
result of low S/N data for A > 25 //m, rather than a physical difference in the environment 
along the hne of sight. 
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4. Calculating the Silicate Extinction Profiles 

In order to compare the laboratory derived silicate minerals to the observed absorption 
spectrum, one needs to calculate extinction profiles from the available optical constants. We 
assume that the size of the silicate grains is much smaller than the wavelength (Rayleigh 
limit). For spherical particles in the Rayleigh limit, the absorption efficiency can be approx- 
imated as 

Q„5, = 4xS5(^j (5) 



where x = 27ra/A and a is the particle size. The cross section per unit volume for absorption 
is 

Cabs 2r> Stt fe-ll 

-^ = 7raQ„,. = -^|-^j (6) 

For a continuous distribution of ellipsoids (CDE), assuming that all shapes are equally prob- 
able, the average absorption cross section per unit volume is 

The above equations are valid for uniform solid particles. To calculate the dielectric function 
for sihcates with inclusions, we used the Bruggeman rule (Bohren & Huffman 1983): 

/ ' ~ +(!-/) ~ - Q (8) 

where eave is the average dielectric constant to be substituted into equations 6 or 7, is 
the dielectric function of the matrix (sihcate) and / is the fiUing factor. For grains will high 
porosity, the filling factor (/) above is replaced with the correction derived by Ossenkopf 
(1991) for spherical inclusions : 

/spheres = /cOS^(7r/). (9) 

The cross-sections per unit volume (Cabs/^) for the silicate minerals and shapes con- 
sidered here are listed in Table 4. 



5. The Silicate Features 
5.1. The Profile of the Interstellar Silicate Features 



In this section we will compare the 9.7 and 18 jim. profiles observed toward the WC stars, 
representative of (local) diffuse ISM dust, with silicate profiles in different environments. 
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5.1.1. The ji Cephei profile 

The emissivity curve for the red supergiant fi Cephei has historically been used to 
represent and model interstellar extinction in the 8-13 /im region (e.g., Roche & Aitken 
1984). The silicates around /x Cep were presumably formed in the outflow of the star. We 
have re-examined the Cep spectrum in light of the identification of photospheric lines in the 
full ISO-SWS spectrum (Tsuji 2000). The photospheric model by Tsuji (2000) shows that 
an absorptive SiO component, whose onset is just shortward of the 9.7 /im silicate emission 
feature, partially underlies it. We first subtract the stellar contimium represented by a 3600 
K blackbody curve (Fig. 3, left panel). The resulting excess emission is fitted with a 250 K 
blackbody curve, representing the underlying dust emission. An emissivity curve is created 
by dividing the fi Cep emission spectrum by the 250 K blackbody curve. The final result 
is shown in the right panel of Fig. 3 along with the absorption spectrum of WR 98a. The 
fjL Cep emissivity curve provides a reasonable representation of the diffuse ISM absorption 
feature. 



5.1.2. The Galactic Center Quintuplet Profile 

The line of sight toward the Galactic center has long been used to probe the properties of 
diffuse ISM dust. Studies of the ISM silicate absorption features toward the central sources 
are comphcated by dust emission (Roche & Aitken 1985; Smith, Aitken & Roche 1990). 
Northeast of the proper center is the Quintuplet Cluster, whose bright sources have also 
been used to study general ISM dust characteristics. The nature of the IR sources remains 
enigmatic, however it is thought that they could be dusty WC stars (Figer etal. 1999; 
Chiar & Tielens 2001; Monieti etal. 2001), much like the sources discussed in this paper. 
Spectropolarimetry of a few of the brightest IR sources (CCS 3-II, III, IV) in the Quintuplet 
Cluster shows that these lines of sight are dominated by pure absorptive polarization across 
the 9.7 /im silicate feature (Smith etal. 2000), making them excellent probes of the ISM 
absorption feature. 

Here we use the Infrared Space Observatory's Short Wavelength Spectrometer spectrum 
of CCS 3 to compare the silicate absorption profiles toward the Quintuplet with those seen 
toward WR 98a. We derived the optical depth spectrum for CCS 3 by fitting a fourth 
order polynomial to the flux spectrum (Fig. 4, top). Fig. 4 (bottom) shows that the relative 
depths of the 9.7 and 18 //m silicate features are similar for GCS 3 and WR 98a. The spectra 
deviate from each other most noticeably in the wavelength region between the two silicate 
absorption features. This could be a result of the difficulty in defining the continuum for 
the WR 98a spectrum longward of the 18 //m feature. Otherwise, the spectra have similar 
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characteristics, including the Tis/rg^j ratio, implying a similar silicate mineralogy for the 
local ISM (WR 98a) and the Galactic center (GCS 3). 

5.1.3. The Asymptotic Giant Branch Star Profile 

Interstellar sihcates are thought to be formed through condensation in the outflows of 
Asymptotic Giant Branch (AGB) stars. As an example of the profiles of sihcates formed in 
such environments, we compare the 9.7 and 18 /im profiles of silicates obtained for diffuse 
ISM with those derived from observations of OH-IR 127.8+0.0 (Fig. 4, bottom). The model 
9.7 /xm extinction profile for the oxygen-rich AGB star OH-IR 127.8+0.0 (Kemper etal. 
2002) is broader on the red wing compared to the ISM (WR 98a and GCS3 I) profiles. This 
additional absorption has been attributed to the presence of crystalline silicates (Demyk 
etal. 2000). The 18 /xm OH-IR star silicate absorption feature peaks at 17.5 /xm (Demyk 
etal. 2000; Kemper etal. 2002), shortward of the diffuse ISM feature which peaks at 18.5 
fjim. Demyk et al. (2000) suggest that the difference in peak wavelength of the 18 //m feature 
is due to the relatively low pyroxene-glass content of the OH-IR star silicates compared to 
those in the ISM (Demyk etal. 2000). 

5.1.4. The Orion profile 

The inner regions of the Orion HII region (Trapezium region) show a 9.7 //m silicate 
emission feature (Forrest, Gillett, & Stein 1975) that is similar to the absorption feature 
observed in the cold dense interstellar medium (e.g., molecular cloud sources; Gillett etal. 
1975). This 9.7 fim profile has been used by Draine & Lee (1984; DL) as constraints on their 
model of interstellar silicates. Figure 5 shows that the DL "astronomical" silicate profile 
does not match the observed central wavelength or width of the 9.7 fira absorption. The 
observed diffuse ISM 9.7 ^ra profile peaks longward of the Trapezium profile. Draine and 
Lee (1984) did not attempt a detailed fit of an observed "20" /xm silicate feature but only 
constrained it to peak at 18//m and to be ~40% of the 9.7 //m feature in strength. The 
feature in the diffuse interstellar medium peaks at slightly longer wavelengths and is slightly 
stronger (deeper) than the DL feature. 
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5.1.5. The Dense Interstellar Medium Profile 

The 9.7 /ill! Trapezium profile has been shown to be more similar to the absorption 
feature observed in dense clouds rather than that in the diffuse ISM. (e.g., Whittet etal. 
1988; Bowey, Adamson & Whittet etal. 1998). Bowey etal. (1998) ascribe the differences in 
the diffuse ISM and dense cloud sihcate absorption profiles to excess absorption, due to a non- 
volatile grain component, on the long- wavelength wing of the dense cloud feature. Similarly, 
Demyk etal. (1999) find that for deeply embedded massive protostars, which probably also 
have dense cloud dust along the line of sight, the silicate absorption profile peaks shortward of 
that in the diffuse ISM and has additional absorption on the long wavelength wing. However, 
for these sources, absorption in the long-wavelength wing is attributed to simple ice species 
(Demyk etal. 1999). This explanation for the profile difference between the diffuse ISM as 
measured toward the WC stars and the Galactic center versus the dense cloud profile would 
not work well for the Trapezium region. First, volatile ices are not expected to survive long 
in the ionized gas of the HII region. In addition, this dust has to be fairly warm in order to 
emit at 10 iim and ice already evaporates at ~ 100 K in the ISM. 



5.2. Comparison with Laboratory Derived Profiles 

The key to identifying the silicate mineral(s) responsible for the observed absorption 
features is finding the combination of mineral(s), shape and porosity that reproduces the 
observed peak wavelengths, widths and relative depths of the observed features. Laboratory 
studies show that the profiles of the silicate features depend on the composition and structure 
of the sihcate material. We will focus here on the results for amorphous sihcates since it has 
been shown that crystalline silicates are not prevalent in the diffuse ISM (Kemper, Vriend 
& Tielens 2004). Silicates are an important component of interstellar dust accounting for 
a major fraction of the interstellar dust volume. Given the observed interstellar silicate 
dust volume, we will focus on silicates based upon the abundant elements, magnesium and 
iron (see also §7). We will only consider here magnesium- iron silicates with an olivine 
or pyroxene stoichiometry. The general chemical formulae for olivines and pyroxenes are 
[(Mg2y,Fe2_2j/)Si04] (with y between and 1) and [Mgj;,Fei_-rSi03], respectively. Generally, 
it is thought that interstellar sihcates are formed in the ejecta of Asymptotic Giant Branch 
stars through reactions with gaseous species (Salpeter 1977). Magnesium silicate compounds 
are calculated to be the first major condensates in a gas with solar system composition cooling 
down from a high temperature (Grossman & Larimer 1974). For the pressures relevant to 
these outflows, at slightly lower temperatures, iron is then expected to react with these 
magnesium silicates rather than condense as separate iron grains, forming iron-magnesium 
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silicates. While the magnesium sihcates are formed at relatively high temperatures and may 
well have a crystalline structure - after all, crystalline forsterite (Mg2Si04) and enstatite 
(MgSiOs) are observed to be abundant in these stellar ejecta, the reaction with gaseous 
iron occurs below the glass temperature of these compounds and, hence, due to kinetics, 
an amorphous structure will result (Tielens & Allamandola 1987; Tielens 1990). While 
strictly speaking the terms olivines and pyroxenes are reserved for crystalline silicates with 
definite mineral structures and hence not appropriate for amorphous silicates, for ease of 
discussion, we will loosely call these amorphous compounds "olivine" and "pyroxene" glasses 
in the subsequent sections. Because of their formation history, these silicates likely have the 
stoichiometry corresponding to olivines and pyroxenes not only on a global level but also on 
a microscopic level and the amorphous nature of these grains reflects small variations in the 
nearest-neighbor bond angle and/or bond distance which destroys all long range order. 

The 10 /im feature in olivine- glass silicates peaks at somewhat longer wavelengths than 
that for pyroxene glasses. At 18 fim, the olivine-glass feature is somewhat sharper and peaks 
at somewhat shorter wavelengths than the pyroxenes. In addition, for both minerals, the 
extinction profile of the 9.7 and 18 /xm features broadens and the strength of the 18 //m 
feature relative to the 9.7 fim feature increases as the Mg:Fe ratio decreases (Dorschner et al. 
1995). The Fe:Mg ratio ranges from to 1.5 and 1 to 1.5 for the pyroxene and olivine glasses, 
respectively, considered by Dorschner et al. (1995). The extinction profiles also depend on the 
shape of the absorbing grains and the porosity of the grains. In the remainder of this section, 
we will make a detailed comparison between the observed silicate profiles and theoretically 
calculated profiles based upon the laboratory measured optical properties (Henning etal. 
1999). 1 



5.2.1. Solid Spheres 

Amorphous olivine-glass is generally used to represent the ISM sihcate absorption fea- 
ture at 9.7 /im. However, while amorphous olivine-glass spheres provide a reasonable match 
to the 9.7 /im absorption feature, the observed 18.5 /im feature is not well matched in either 
central wavelength or relative depth. As shown in Fig. 6 (top left panel), the spectrum of 
olivine-glass spheres peaks near 17 /im and is too deep relative to the observed feature. The 
spectra of amorphous pyroxene-glass spheres peak near 9.3 /xm and therefore do not pro- 



^ Optical constants for amorphous olivine and pyroxene glasses from the Laboratory Astrophysics 
Group of the Astrophysical Institute and University Observatory in Jena (http://www.astro.uni- 
jena.de/Laboratory/Database/odata.html) 
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vide a good match to the observed 9.7 //m feature. On the other hand, the relative depths 
of the 9.7 and 18 features, as well as the central wavelength of the 18 /im feature are 
well- reproduced by amorphous pyroxene glasses (Fig. 6, top left panel). 

Thermodynamic equilibrium calculations show that, for temperatures less than about 
1300 K^, solid forsterite is converted into solid enstatite with excess gaseous SiO (Grossman 
and Larimer 1974). The extent to which this reaction proceeds before grain formation freezes 
out in the stellar ejecta depends on the timescale and density/temperature structure in the 
outflow and, hence, is difficult to constrain observationally. For this reason, we attempted to 
reproduce the astronomical features by a mix of these two components (assuming separate 
populations). The bottom left panel in Fig. 6 shows that the match using an arbitrary 
pyroxene- glass/olivine-glass mixture is improved over considering only one silicate species. 
However, the match between the pyroxene- glass/olivine-glass "mixture" and the observed 
spectrum is less than ideal between 11 and 15 /im. In addition, the peak wavelengths of 
both features are somewhat blueward of the observed feature. One way to shift the peak 
wavelengths redward is to incorporate vacuum into the grains. We discuss such porous grains 
below. 



5.2.2. Porous Spheres 

Observations of dust in dense clouds suggest that grain growth is dominated by coagu- 
lation (Jura 1980), thus grains are likely to be porous. Porous grains have larger extinction 
cross sections relative to their solid counterparts and therefore require less of a given element 
to account for, e.g., the silicate interstellar absorption features. Mathis (1998) argues that 
the cosmic abundance requirements can be met if 25% of the grain volume is vacuum. An 
upper limit of 60% on the porosity is given by the width of the interstellar polarization law 
(Wolff, Clayton, & Meade 1993). 

We used the Bruggcman rule (Bohrcn & Huffman 1983) to calculate the extinction 
profiles for porous silicates, assuming that the vacuum is in the form of spherical inclusions 
(§4). In general, increasing the vacuum volume has the effect of decreasing the 9.7 fim to 
18 /im ratio and shifting the peak absorption of both bands redward. Fig. 6 (top right 
panel) demonstrates how the extinction profiles of pyroxene-glass and olivine-glass spheres 
are modified by 50% vacuum, pyroxene-glass spheres that are 50% porous do a fair job at 
reproducing the astronomical spectrum, although they provide too much extinction on the 
blue side of the 9.7 /im feature. Porous olivine-glass spheres (50% vacuum) do not reproduce 



^The exact temperature where this conversion sets in depends somewhat on the pressure. 
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the observed spectrum, the biggest problem is the 18 //m to 9.7 //m band ratio that is too 
high to match the observations. The extinction profile of a mixture of porous (50% vacuum) 
olivine-glass and pyroxene-glass spheres is shown in the bottom right panel of Fig. 6. In this 
case, the 9.7 fim band is somewhat broader than the observed feature and the 18.5 /xm to 
9.7 iim band ratio is overestimated, although not significantly. 

5.2.3. Continuous Distribution of Solid Ellipsoids 

A continuous distribution of ellipsoids (CDE) where each shape occurs with equal prob- 
ability could reflect the characteristic grain population in the ISM. Figure 7 (left panels) 
shows the calculated extinction profiles, using equation 7, for solid CDE. The extinction 
profiles for CDE arc somewhat broader and are shifted rcdward relative to those for spheres. 
As shown in Fig. 7, the extinction profile for olivine-glass CDE is broader on the red side of 
the 9.7 /im profile, and produces a stronger 18.5 /im feature, relative to the observed pro- 
file. On the other hand, the profiles calculated for pyroxene-glass CDE match the relative 
strengths of the 9.7 and 18.5 /xm features well, but are broader than the observed 9.7 iim 
profile on the blue side. Co-adding the extinction profiles of pyroxene-glass and olivine-glass 
CDE improves the match to the wings of the 9.7 fim profile, however, the relative strength 
of the 9.7 /im feature is underestimated (Fig. 7, bottom left panel). 

5.2.4- Continuous Distribution of Porous Ellipsoids 

We use the Bruggeman rule again to calculate the extinction profile for porous CDE 
with 50% porosity. The extinction profiles for porous CDE result in a somewhat increased 
18 /xm to 9.7 fj,m feature ratio compared to solid CDE. In addition, both profiles are broader 
than those for solid CDE. 



5.2.5. Alumino silicates 

Aluminum is heavily depleted in the interstellar medium (Whittet 2003), and as such, 
is likely to be incorporated into silicate dust. Aluminosilicatcs were considered for circum- 
stellar silicates around oxygen-rich stars because they provide much-needed opacity in the 
trough region between the 9.7 and 18 ^m. features (Mutschke et al. 1998). We considered Fe- 
and non-Fe- containing aluminosilicatcs for the diffuse ISM absorption features. Extinction 
profiles for porous and solid spheres and CDE were compared with the observations based 
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on the optical constants by Mutschke etal. (1998). Fitting the observed sihcate profiles 
with aluminosilicates (with and without Fe) with Si:Al equal to 2:1 requires more Al than 
available in the ISM independent of grain shape and porosity. The Fe-containing aluminosil- 
icates provide more opacity in the trough region compared to aluminosilicates with no Fe. 
In general, for the alumino-silicates considered, the relative depths of the observed absorp- 
tion features are matched by the extinction profiles for porous spheres, or solid or porous 
CDE; sohd spheres result in a larger 9.7/im/18 fim ratio than observed. The observed 18 
/im profile shape is also well-matched by the extinction profiles for porous spheres, or solid 
or porous CDE. Aluminosilicates look less promising in the spectral region of the observed 
9.7 /im absorption feature. Here, the calculated extinction profiles do not provide a good 
match to the observed peak wavelength and the profile width, simultaneously. Thus, it is 
unlikely that aluminosilicates can account for the observed silicate profiles on their own, but 
they may still be present as a component of the interstellar silicate dust. 



5.2.6. Iron Oxides 



Theoretical calculations under conditions appropriate for the cooling solar nebula (e.g., 
high pressures as compared to stellar ejecta) predict that iron oxide forms by oxidation of 
metallic iron with excess H2O at about 600 K (Barshay & Lewis 1976). The fate of con- 
densing iron is less clear in the cjccta of asymptotic red giant branch stars. Depending on 
temperature, pressure, and oxygen partial pressure, iron condensation can shift from metallic 
iron, to iron- magnesium silicates, to ferrous iron (FeO, wuestite), mixed ferrous- ferric oxide 
(Fe304; magnetite), to ferric iron (Fe203; hematite) in thermodynamic equilibrium calcu- 
lations. Of course, kinetic factors may hamper the formation of some of these compounds 
(eg., metallic iron). In interstellar and circumstellar medium settings, the simplest iron ox- 
ide, ferrous oxide (FeO), has at times been considered as an important contributor to the 
extinction in the infrared (Hennning et al. 1995; Demyk ct al. 2000; Kemper et al. 2002). The 
presence of circumstellar magnesio- wuestite grains has been inferred from the 19 /im emission 
feature in the spectra of AGB stars with low mass loss rates (Cami 2002; Posch etal. 2002; 
Heras & Hony 2005). However, those grains are considered to be largely magnesium- rich. 
On Earth, wuestite (FeO) is considered an important constituent of the lower mantle and 
natural samples of wuestite are available. In the laboratory, wuestite is routinely made using 
a gel method and its optical properties have been measured for astronomical purposes. Of 
course, in the presence of excess oxygen, ferrous oxide burns explosively at temperatures 
exceeding 600 K, but - in view of kinetic considerations - that may not be very relevant 
for interstellar conditions (where temperatures seldomly reach such high values), anyway. 
Here, in view of these discussions in the literature, we merely examine the effect of small 
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FeO inclusions in silicate grains on the 9.7 and 18 //m silicate profiles if they were present 
in the ISM. 

Spherical inclusions of FeO are considered here for interstellar silicates using the Brugge- 
man rule. As discussed by Henning etal. (1995), the presence of FeO in sihcates has the 
affect of increasing the 18//m/9.7//m ratio and broadening both bands. The 9.7 /im and 18 

/xm features shift blueward and redward, respectively, as the volume of FeO inclusions is 
increased. Thus, adding FeO (spherical) inclusions to solid olivine-glass or pyroxene-glass 
spheres would worsen the match to the interstellar spectrum. For solid olivinc-glass spheres, 
FeO inclusions increases the depth of the 18 fira band, and for solid pyroxene-glass spheres, 
the 10 yUm band is shifted to even shorter wavelengths. Similarly, adding FeO to porous 
spheres does not improve the fit. Therefore, iron oxides are probably not a significant grain 
component in ISM dust. 

5.2.7. Summary 

Considering the parameter space investigated here, we find two equally good matches to 
the 9.7 and 18 /xm silicate absorption profiles. Our criteria for a good match are as follows: 
the extinction profile should fit the relative depths, the central wavelengths, and the wings 
of both the 9.7 and 18 iim absorption features. One good match is achieved by co-adding the 
extinction profiles of olivine-glass and pyroxene-glass spheres with 50% porosity. An equally 
good match is given by the co-addition of pyroxene-glass and olivine-glass solid CDE. Both 
these fits require a greater contribution by mass from pyroxene-glass relative to olivine-glass. 
This differs from the conclusion reached by Kemper, Vriend & Tielens (2004) that olivine 
glasses account for the majority of the mass of diffuse ISM silicates. However, these authors 
consider only the 9.7 /im feature which we also find is best matched with olivine glasses if 
considered by itself. 

Several of the calculated extinction profiles shown result in unsatisfactory fits to the 
observed silicate profiles. While the extinction profile for solid olivine-glass spheres matches 
the observed 9.7 /^.m silicate profile well, it does not match the central wavelength or relative 
depth of the 18 /im feature (Fig. 6, top left). On the other hand, solid pyroxene-glass spheres 
match the 18 /im feature, but not the 9.7 /im feature. Co-adding the pyroxene-glass and 
olivine-glass profiles for sohd spheres does not improve the fit (Fig. 6, bottom left). Fig. 6 
(top right) shows that adding porosity improves the fit for pyroxene-glass spheres, though 
the blue wing of the 9.7 /xm feature is not well-matched. Porous olivine-glass spheres do not 
match either the 9.7 or 18 /im profiles. Neither the solid CDE or porous CDE extinction 
profiles for pyroxene or olivine glasses on their own result in satisfactory fits to the observed 



-14- 



silicate features (Fig. 7). 

Aluminosilicates are good candidates for the opacity in the trough region between the 
two sihcate absorption features. Many of the calculated extinction profiles also match the 
18 /im profile and the relative depths of the features well. However, these minerals do not 
provide a good match to the 9.7 /im profile shape or peak wavelength. Thus, we conclude 
that these minerals arc unlikely to account for the observed features on their own, but could 
partially contribute to the observed absorption. We also find that the presence of iron oxides 
cannot be accounted for in the observed ISM spectrum. 



6. The Lifecycle of Interstell£tr Silicates 

The hfecycle of interstellar sihcates starts with the nucleation and chemical growth at 
high densities and temperatures in the ejecta from stars such as Asymptotic Giant Branch 
stars and supergiants. This ejected material is rapidly mixed with other gas and dust in 
the interstellar medium. In the interstellar medium, dust cycles many times between the 
intercloud and cloud phases until it either is destroyed by fast (~ 100 km/s) supernova 
shocks or is incorporated into newly formed stars or planetary systems. 

Infrared observations with the Infrared Space Observatory have revealed that crystalline 
silicates are abundant (~ 15% by volume) in the initial stages - the ejecta from stars - as 
well as in the last stages of its life cycle - circumstellar disks around Herbig AeBe stars and 
T-Tauri stars - but are absent in the middle portion of this evolutionary scenario (Sylvester 
etal. 1999; Kemper etal. 2004, 2005; Malfait etal. 1999; Waters & Waelkens 1998). Indeed, 
the smooth profile of the 9.7 /im silicate feature implies that the crystalline-to-amorphous 
sihcate fraction is less than 1% in the interstellar medium (Kemper et al. 2004, 2005). Likely, 
the crystalline silicates injected in to the interstellar medium by stars are amorphized by high 
energy galactic cosmic rays on a timescale (70 million years; Bringa et al. 2005) which is short 
compared to interstellar dust evolution timescales (~ 500 — 4000 million years; Jones et al. 
1996). The observed high abundance of crystalline silicates in circumstellar disks surrounding 
young stars must then reflect subsequent processing of grains in these environments. Recent 
mid-infrared spectroscopic interferometry on AU-scale-sizes has revealed the presence of a 
strong gradient in the crystallinity of silicates in the disks of the three objects investigated 
(van Boekel etal. 2004). Likely this reflects the rapid annealing in the hot inner regions of 
these disks coupled with turbulent diffusion outward of the crystalline grains. 
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7. Elemental Abundances 

The fraction of the elements locked up in interstellar dust grains provides a key test 
for any realistic dust model (e.g. Snow & Witt 1996; Mathis 1996). Interstellar abundances 
have recently been reviewed by Sofia & Meyer (2001) who show that young (< 2 Gyr) F and 
G stars and the Sun provide the best representations of the interstellar abundances. 

We calculate the silicate abundance per unit hydrogen atom for both solid and 50% 
porosity silicate (olivine-glass and pyroxene-glass) spheres and CDE, and tabulate the results 
in Table 5. The final row in Table 5 hsts the silicate abundance using a mixed population 
of olivine-glass and pyroxene-glass grains. We assume an average total hydrogen column 
density per unit extinction (< N-a/Av >) of 1.9 x 10^-*^ cm~^ mag~^, with a typical scatter of 
less than 30% over a wide range oi E{B - V) (Bohhn et al. 1978) . We take Ay /T9.7 = 18 ± 1 , 
determined for the interstellar medium (Roche & Aitken 1984). Then, the abundance of 
silicate atoms relative to hydrogen locked up in solid form is, 

-^silicate 1 /^silicate 1 /' 1 n^ 

A^H ~ Cabs/V 1.66 X 10-24M,aicate 18.5 X 1.9 X 1021 ^ > 

where PsiUcate a-nd Mgiiicate a-re the specific density and the atomic mass per silicon atom 
of the silicate material. Our discussion on elemental abundances will be based on the evalua- 
tion by Sofia & Meyer (2001). They conclude that the atomic Si available to be incorporated 
into dust is between 3.44 — 3.99 x 10^^ per H atom, enough to alleviate the Si underabun- 
dance problem proposed by Snow & Witt (1996), for many (but not all) silicate dust models. 
Using solid olivine-glass and/or pyroxene- glass spheres to account for the observed silicate 
profiles requires 5 to 49% more than the available Si. Similarly, modeling the silicate ab- 
sorption features with solid pyroxene-glass CDE requires somewhat more than the available 
Si. Assuming that interstellar silicate grains are porous allows the interstellar features to be 
accounted for without using up the available Si. There is enough available Mg and O to be 
incorporated into silicate grains. The elemental abundance of Fe, on the other hand, poses 
a problem for olivine-glass grains unless porous grains are used (Table 5) . 

Comparison of elemental abundances in the cold neutral medium as compared to the 
warm neutral medium suggests that a fraction of 0.3 of the elemental Si is locked up as a dust 
component which is not as refractory as silicate materials and is more readily destroyed by 
interstellar shocks (Sembach & Savage 1996; Tielens 1997). If we assume that we only have 
2.4 — 2.8 X 10^^ Si per H-atom available for interstellar silicate dust, all models would fail 
the test (Table 5). However, porous spheres and CDE models would come the closest. The 
small discrepancy may just reflect a small difference in the elemental abundance of silicon 
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in the direction of the inner Galaxy due to the elemental abundance gradient of the galaxy. 
From other elements, we estimate that this may amount to an increase in the Si abundance 
of a factor of 1.2 over the typical distance (~2 kpc) for which the local value of Ay/rg.y was 
determined. 



8. Other dust components 

Silicon carbide is produced in the outflows of carbon-rich low-mass asymptotic giant 

branch stars (Blanco ctal. 1994; Groenewegen 1995; Speck ctal. 1997) and is incorporated 
into meteoritic dust (Anders & Zinner 1993 and references therein). It should be prevalent 
in ISM dust, but limits of < 5% have been placed on its abundance (Whittet, Duley, & 
Martin 1990). 

There is no evidence for SiC grains in the observed spectra. We put a limit on the 
amount of SiC present in the ISM using the (uncorrected) extinction efficiencies given by 
Borghesi et al. (1987) for /3-SiC since it peaks near the observed carbon-star emission feature. 
Its extinction efficiency is ((5cxt/a)max = 3.32 x 10*^ cm~^. Assuming a specific density s = 3.2 
g cm^'^, the mass absorption coefficient, k, is calculated to be 7780 cm^ g^^ To calculate 
the column density of Si atoms contained in SiC grains, we use the equation (Whittet et al. 
1990) 



iv^(Si) = -Al^ (11) 

where fx is equal to 0.7, the mass-fraction of Si in SiC. The optical depth is calculated 
from the spectra of WR 98a and WR 112. We do not use the WR 104 spectrum since, 
at wavelengths longer than 11 /im, it significantly deviates from the mean of the other 3 
spectra. For WR 118, olivine-glass and pyroxene-glass already slightly overcompensate for 
the absorption around 11.3 /im. 

We use a x^-fitting routine to determine the maximum allowable depth of /3-SiG in 
the WR 98a and WR 112 spectra using a three-component fit (/3-SiG, olivine-glass and 
pyroxene-glass solid spheres). We determine Th.s < 0.016 and 0.019 for WR 98a and WR 
112, respectively. Thus, the limit on the fraction of Si in SiC with that in silicates is 0.03 and 
0.04 for WR 98a and WR 112, respectively. The results of our calculations are tabulated in 
Table 6. 
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9. Extinction Curve 

The extinction curve in the infrared is generally accepted to be uniform between 0.9 and 

5 nm. Based on an extensive compilation of photometric data for dense clouds and diffuse 
ISM sightlines, Martin & Whittet (1990) suggest a universal interstellar extinction curve 
from 0.35 to 5 /xm, represented by a power law, Ax/ Ay — X~^^^. Recently, Indebetouw et al. 
(2005) used photometric data from 1.25 to 8 //m of the Galactic plane from both the Spitzer- 
GLIMPSE program and 2MASS to investigate the interstellar extinction curve and extend 
its validity to 8 fim. They deduce the extinction law, logfAA/A^] = 0.61 — 2.221og(A) + 
1.21[log(A)]^, that is valid for the dense and diffuse ISM within the solar circle. These authors 
show that their extinction curve provides a reasonable match to the Lutz (1999) Galactic 
center (GC) values that were deduced from hydrogen recombination lines, demonstrating 
that, at least in the near-IR, the extinction in the local ISM and the GC are similar. We 
plot both data sets in Fig. 8, renormahzing the Lutz (1999) points to the K-band extinction, 
assuming Ay — 29 magnitudes and A^ — 3.28 magnitudes (Figer, McLean & Morris 1999). 
Given the good agreement of the Indebetouw etal. (2005) and the Lutz (1999) extinction 
points, we combined them and carried out a least-squares fit assuming the same functional 
form as Indebetouw etal. (2005) and find log[Ax/AK] = 0.65 - 2.401og(A) + 1.34[log(A)]2. 
We use this equation to represent the continuum extinction, assuming that at A > 8 /im, the 
continuum remains constant and the silicate extinction dominates. 

It is well-known that there is more silicate per unit of visual extinction toward the GC 
compared to the local ISM (e.g., Roche & Aitken 1984, 1985). Thus, longward of 8 /im, 
we need to construct individual extinction curves for the local ISM and the GC. To create 
the "GC extinction curve" we first normalize the silicate optical depth spectrum to the 
average visual extinction to silicate optical depth ratio, Ay/rg.y = 9, determined by Roche 

6 Aitken (1985) for the Galactic Center. The spectrum was then normalized to Ak using 
the Figer etal. (1999) Ay and Ak values for the Galactic Center (see above). Finally, this 
curve was added to our continuum extinction curve assuming that the "continuum" remains 
constant longward of 8 //m. The final GC extinction curve is shown in Fig. 8 (left panel) and 
partially tabulated in Table 7. A complete downloadable version of the table is available in 
the electronic edition of the Journal. 

To create a complete extinction curve for the local ISM, we carry out the same set of 
steps as for the GC, but this time we use Ay/rg./ = 18, appropriate for the local ISM (Roche 
& Aitken 1984; Rieke & Lebofsky 1985). In order to normalize this to the extinction in 
the K band, we use Ak/Av = 0.09 (Whittet 2004). Finally, we add the normalized WR 
98a silicate profile to the continuum extinction. The complete local interstellar extinction 
curve is shown in Fig. 8 (right panel). A partial list of tabulated values is shown in Table 7; 
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a complete downloadable version of the table is available in the electronic edition of the 
Journal. 



10. Summary and Discussion 

In this paper, we examined the silicate absorption features observed in the diffuse ISM 
toward WC-type WR stars. As earher work has shown, the observed 9.7 //m feature is most 
similar to the /i Cephei profile, and differs from the Trapezium profile. The latter is used 
by Draine & Lee (1984) to compute their "astronomical" silicate profile, and is found to be 
more similar to the silicate absorption features observed in dense clouds. We have produced 
infrared extinction curves appropriate for the GC and local ISM. These two curves, which 
range from 1.25 to about 30 /xm are available for download in the electronic version of the 
Journal. 

We find that pyroxene-glass silicates dominate by mass in the diffuse ISM silicates. 
Pyroxene-glass silicates are also found to dominate the spectra of YSOs (Dcmyk etal. 1999; 
Dorschner etal. 1995 and references therein; Jaeger etal. 1994). In contrast, for evolved 
stars, it has been estimated that no more than 10% of the amorphous silicates (by mass) are 
pyroxene-glass. It is unlikely that this reflects a genuine conversion of olivine to pyroxene 
glasses at the low temperatures of the diffuse interstellar medium. Rather, we attribute this 
difference to an incomplete sampling of the stellar sources that inject silicate dust into the 
interstellar medium. In particular, supernovae are conceivably an important source of inter- 
stellar dust. However, because of their scarcity in the solar neighborhood, infrared spectra 
of SN dust ejecta are virtual non-existing. With the launch of Spitzer and, particularly, its 
sensitive infrared spectrometer, IRS, this situation will hopefully quickly be remedied. 

We note that, in our fits, interstellar amorphous silicates contain an appreciable frac- 
tion of iron. In contrast, crystalline silicates injected by late-type stars are in the form of 
forstcritc and enstatite - the magnesium rich end members of the olivine (Mg2Si04) and 
pyroxene (MgSiOa) families. Theoretically, the condensation behavior of iron in the ejecta 
of late-type stars is somewhat different than for the Solar Nebula. Indeed, at the low pres- 
sures of AGB winds, iron may not condense out until after the magnesium rich silicates have 
formed. In that case, because of the high super-saturation required for nucleation of metallic 
iron, rather than condensing as separate iron-metal grains, iron may react with these sili- 
cates to form mixed iron-magnesium silicates. Because this reaction occurs below the glass 
temperature of these silicates, this reaction would be expected to lead to the formation of 
glassy silicates (Tielens 1989). Within this scheme, the crystalline silicates would form at 
higher temperatures, well above the glass temperatures, by direct condensation from the gas 
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phase. 

As remarked above, the profile of the 9.7 iim feature in dense clouds and toward YSO's 
has additional absorption in the wing compared to that of the diffuse ISM. While the profiles 
differ, the amorphous silicates in these different environments are generally assumed to be 
miner alogically similar. The additional absorption in dense clouds is generally attributed to 
a non-volatile dust species, whereas for the YSO's the additional absorption is accounted for 
by simple ice species. Such fits are however not unique. Some processing of silicates may 
occur in the protoplanetary disks around protostars, resulting in a different mineralogy. This 
material may then be deposited into the surrounding molecular cloud through the powerful 
protostellar winds and hence the general interstellar medium (Tielens 2003). However, how 
much processing occurs in these environments and how well this material is mixed into 
molecular clouds and the general interstellar medium is presently not well understood. 

Recently, genuine silicate Stardust grains have been isolated in interplanetary dust par- 
ticles(IDPs) (Messenger etal. 2003). Some of these grains are crystaUine while others are 
amorphous. While most presolar grains found in this study of IDPs contain abundant O, Si, 
and Mg, their mineralogy is currently unknown. Only one grain is identified as forsterite. 
Moreover, some interstellar dust may have a very mundane isotopic (eg.. Solar abundances) 
composition and these grains would be difficult to separate from Solar system condensates 
(Tielens 2003). Nevertheless, such studies hold the promise of direct information on various 
relevant aspects of interstellar silicates. We note that GEMS, Glass with Embedded Metals 
and Sulphides, arc an abundant component of anhydrous IDPs (Bradley 1994). The de- 
tailed morphology and structure of these GEMS betrays an extensive ion-irradiation history 
(Bradley 1994). While this may reflect long exposure to the ISM, this irradiation may also 
have occurred in the early solar system environment (Tielens 2003; Westphal and Bradley 
2004). The 9.7 //m band in GEMS-rich thin sections of IDPs is broad and featureless and, in 
that respect, resembles the 9.7 /xm features observed towards YSOs, the Trapezium cluster, 
and fi Cep (Bradley etal. 1999). Of course, as demonstrated in this paper, these interstellar 
9.7 fim features differ in detail. Nevertheless, these initial laboratory studies of solar system 
materials, which may well be directly related to interstellar silicate grains, are very promis- 
ing. It is clear that many questions remain on the detailed structure and composition of 
interstellar sihcates. 
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Table 1. Observing parameters 



Source 


TDT 


AOT 


WR 98A 


09401206 


AOTOl 


WR 104 


09901207 


AOTOl 


WR 112 


10201908 


AOTOl 


WR 118 


10802509 


AOTOl 



Table 2. Model continuum parameters for the WR stars 



Source 


T 


Ti 




K 


K 


WR 98A 


790 


180 


WR 104 


840 


110 


WR 112 


680 


90 


WR 118 


1020 


150 



Table 3. Silicate Optical Depths and Visual Extinctions 



Source 


T9.7 


^18.5 






WR 98A 


0.78 


0.39 


2.0 


12.42 


WR 104 


0.38 


0.28 


1.4 


6.50 


WR 112 


0.63 


0.32 


2.0 


11.03 


WR 118 


0.78 


0.45 


1.7 


11.20 



^A(V) from van der Hucht 2001. 
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Table 4. Calculated Silicate Band Strengths (Cext/1^ [10^ cm-^]) 



Amorphous Silicate 


Solid Spheres 


Solid CDE 


Porous Spheres^ 


Porous CDE' 


OUvine (MgFeSi04) 


0.880 


0.926 


0.616 


0.620 


Pyroxene ((Mgo.sFeo.sSiOs) 


0.983 


0.997 


0.666 


0.656 



^Fifty percent porosity. 



Table 5. Sihcon Abundance Relative to Hydrogen (xlO ^) 



Amorphous Sihcate Solid Spheres Solid CDE Porous Spheres^ Porous CDE^ 



Ohvine (MgFeSi04) 


4.20 


3.99 


3.00 


2.98 


Pyroxene (Mgo.sFeo.sSiOa) 


4.80 


4.74 


3.55 


3.60 


Olivine + Pyroxene 


4.34 


3.89 


3.14 


3.07 




(3.28+1.06) 


(1.52+2.37) 


(0.87+2.27) 


(0.27+2.80) 



^Fifty percent porosity. 
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Table 6. Limit on SiC in ISM 



Source 


Olivine Pyroxene 


N(total silicate) 


N(/3-SiC) 


N(Sisic)/N(Sisi,icates) 




xio-Vh 


xlO^^ cm 


-2 




WR 98a 


3.02 1.15 


98.4 


3.07 


0.03 


WR 112 


2.86 1.34 


88.0 


3.65 


0.04 



Table 7. Extinction curve for local ISM and Galactic Center^ ^ 



Wavelength 
/im 


Ax/Ak 
local ISM 


Ax/Ak 
Galactic Center 


1.240 


2.724 


2.724 


1.340 


2.315 


2.315 


1.440 


2.002 


2.002 


1.540 


1.759 


1.759 


1.640 


1.565 


1.565 


1.740 


1.408 


1.408 


1.840 


1.279 


1.279 


1.940 


1.171 


1.171 


2.040 


1.081 


1.081 


2.140 


1.004 


1.004 



^The complete version of this table is in the 
electronic edition of the Journal. The printed 
edition contains only a sample. 

''The extinction curves for the local ISM and 
GC are equal at wavelengths less than 8 /im. 
Longward of 8 /im, the extinction curves for 
these environments differ; see §4 and Fig. 3. 
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Wavelength (^im) 

Fig. 1.— 2.5-45 fim ISO-SWS AOTl spectra for heavily extincted WR stars WR 112, WR 
118, WR 104, and WR 98A. All spectra show the 9.7 and 18 fim silicate absorption. Also 
shown are photometry points from Williams etal. 1987 (where available), normalized to the 
ISO spectrum in the L band. The continuum model is described in the text. 
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Fig. 2. — Optical depth absorption spectra resulting from model continuum fits shown in 
Figure 1. The order of spectra from top to bottom is WR 104, WR 112, WR 118 and WR 
98a in both panels. [Top] The 7 to 40 iim region. The spectra are offset from each other 
for clarity. [Bottom] The 7.5 to 13 //m region of the sihcate feature. The spectra have been 
normalized to rg.y = 1.0 and each spectrum is offset by 0.2 optical depth units from the 
spectrum below. 
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Fig. 3.— [left] The dereddened {Ay = 1.5) ISO-SWS spectrum of fi Cep shown with a 3600 
K blackbody representing the stellar continuum, [right, dashed line] The n Cep emissivity 
curve derived from the spectrum at the left assuming = 250 K for the underlying dust 
emission, compared with the silicate absorption feature of WR 98a [solid line]. 
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Fig. 4. — [Top] ISO-SWS spectrum of GCS 3-1 with fourth order polynomial continuum used 
to derive the optical depth spectrum in the panel below. [Bottom] Extinction model profile 
for OH-IR 127.8+0.0 (dashed line; Kemper et al. 2002). Extinction profile for the Galactic 
Center Quintuplet source GCS 3-1 (dotted line) . Extinction profile for WR 98a (solid line) . 
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Fig. 5. — The extinction profiles for Drain & Lee "astronomical silicates" compared to the 
WR 98a absorption profile. 
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Fig. 6. — The absorption profile for WR 98a normalized to unity at 9.7 fim compared 
with the extinction profiles for amorphous olivine-glass (dotted line) and pyroxene-glass 
(dashed line) spheres. Solid [left panels] and porous [right panels] spheres are considered. 
The porous spheres contain 50% vacuum. The bottom panels show the effect of combining 
separate populations of olivine-glass and pyroxene-glass spherical grains; the solid smooth 
line represents the total absorption profile. 
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Fig. 7. — The absorption profile for WR 98a normalized to unity at 9.7 iim compared with 
the extinction profiles for amorphous olivine-glass (dotted line) and pyroxene-glass (dashed 
line) CDE. Sohd [left panels] and porous [right panels] CDE are considered. The porous CDE 
contain 50% vacuum. The bottom panels show the effect of combining separate populations 
of olivine-glass and pyroxene-glass CDE; the solid smooth line represents the total absorption 
profile. 
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Fig. 8. — Interstellar extinction curves for the local ISM (left) and the Galactic Center 
(right). Circles represent measured extinction for the Galactic plane (Indebetouw etal. 
2005. Triangles represent the extinction measured toward the Galactic Center (Lutz 1999). 
The smooth line from 1.25 to 8 fim is a least-squares fit to the measured extinction points. 
The silicate curves are representations of the silicate features for the local ISM (WR 98a; 
left) and the Galactic Center (GCS3; right). 



